Phase Transition in a Self-repairing Random Network 



m 
o 
o 

(N 

c 



I 

CZ3 



o3 



C 

O 

o 



> 

O 

CO 

o 

-I— > 
03 



T3 

o 
o 



X 



A. S. Ioselevich, D. S. Lyubshin 

Landau Institute for Theoretical Physics, Russian Academy of Sciences, Kosygina str.2, 117940 Moscow, Russia 

(Dated: February 2, 2008) 

We consider a network, bonds of which are being sequentially removed; that is done at random, 
but conditioned on the system remaining connected (Self- Repairing Bond Percolation SRBP). This 
model is the simplest representative of a class of random systems for which forming of isolated 
clusters is forbidden. It qualitatively describes the process of fabrication of artificial porous materials 
and degradation of strained polymers. We find a phase transition at a finite concentration of bonds 
p = p c , at which the backbone of the system vanishes; for all p < p c the network is a dense fractal. 

PACS numbers: 



Properties of networks, bonds (or nodes) of which are 
being removed randomly, have been intensively studied 
during past 50 years. In the standard formulation of 
the problem the bonds (or sites) are removed totally at 
random, and the percolation transition takes place at a 
certain concentration p c of remaining bonds (or concen- 
tration x c of remaining sites). For p > p c there exists 
an "infinite" cluster that contains a finite fraction of all 
bonds in the system and spans homogeneously through 
the entire network. At p < p c the infinite cluster does not 
exist and only finite ones are present. The percolation 
phase transition and corresponding critical phenomena 
are well studied (see, e.g., 

However, there are many physical systems which can- 
not be described by the standard percolation theory. In 
particular, there are important cases when finite clus- 
ters can not appear at all. As an example, consider a 
technologically important process of pore forming (i.e., 
fabrication of a porous material; see, e.g., 0)- It can 
be viewed as gradual removal of grains of a pore-former 
(carbon, which can be burned out, or a soluble polymer) 
from a mixture of the pore-former with grains of a ma- 
trix material (a metal). Due to mechanical instability 
of finite clusters, they immediately fall down onto the 
surrounding matrix and stick to it. Thus at any stage 
of the process the remaining grains form a single "infi- 
nite cluster" and the percolation transition is impossible. 
The properties of this single cluster are nontrivial: it ap- 
pears that there is a topological phase transition at a 
finite concentration of remaining grains x c , below which 
the system becomes "cracked" and its mechanical and 
conducting properties degrade catastrophically. 

So far we were not able to demonstrate the existence 
of the above phase transition by analytic methods, our 
qualitative arguments and results of numerical studies 
will be published elsewhere Q. In the present paper we 
introduce a simplified model, which allows rigorous anal- 
ysis. The model is very similar to the standard bond 
percolation: starting from a full lattice, at each step one 
of the remaining bonds is randomly chosen for removal. 
But after its removal the system is checked for the ex- 
istence of finite clusters: if such ones are present, then 



the removed bond is restored (i.e., the last removal is 
cancelled) and the process goes on to the next step. It 
seems natural to call this model a "self-repairing bond 
percolation" (SRBP). 

Apart from being relevant to pore-forming, the SRBP 
model may also be viewed as a model for polymer degra- 
dation (see, e.g., Consider a random network con- 
sisting of irregularly cross- linked polymer chains. Sup- 
pose that this system is subjected to random external 
perturbation (e.g., UV-radiation) that can destroy the 
cross-links. The radiation damage may be repairable: 
attraction between individual chains tends to reestablish 
the damaged link. However, sometimes that appears to 
be impossible, since internal strains in the chains may 
drive the two chains apart as soon as the link between 
them is damaged. Thus it seems reasonable to assume 
that all strained links are vulnerable to radiational dam- 
age, while unstrained ones are "immune" to it. Of course, 
finding out which links in a random network are strained 
and which are not is a formidable task. But in any case 
the links which are the only bridges connecting otherwise 
isolated clusters are never strained. These links will be 
repaired after possible removal, in accordance with the 
definition of the SRBP model. 

Let the fraction of the remaining bonds be p. We will 
be interested in average properties of the system as func- 
tion of p. In particular, we will study the conductivity 
a(p) and "the minimal chemical path" £(R,p), the latter 
being an ensemble-averaged length of the shortest path 
going via bonds and connecting two points separated by 
euclidean distance R. 

The first obvious observation about the SRBP model 
is that there exists some minimal possible p = puce = f 
[z being the coordination number of the lattice) at which 
the process of bond removal stops: for a connected graph 
one necessarily has p > ptroc- At p — pt rcc the re- 
maining bonds constitute a spanning tree (ST), a con- 
nected graph with no cycles and all lattice sites as ver- 
tices. As it is well known, the probability of generat- 
ing a given ST at the end of our process is related to 
the Minimal Spanning Tree (MST) problem (see, e.g., 
6J). In particular, the minimal chemical path is fractal: 
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FIG. 1: A typical configuration of black ang grey bonds at 
the critical concentration for 2D square lattice. 
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£(R,P = Ptrecj 

for a tree one has a(p = ptree) 

We will show that actually the minimal chemical path 
is fractal and the specific conductivity of the system is 
zero not only at p = pt iee , but also within a finite interval 
Ptrcc < P < Pc ■ The corresponding phase we will call the 
"tree-like phase" , in contrast to the "solid phase" existing 
at p c < p < 1 . 

To prove the above statement, we use a mapping to 
the standard percolation. Suppose that initially all the 
bonds of the system are black. If at some step a removal 
of a certain bond must be cancelled, we restore the bond 
but change its color to grey. Then for any fraction p 
of remaining bonds we have fractions b = b(p) of black 
and g = p — b{p) of grey bonds remaining, where b(p) 
is a certain monotonically increasing function with the 
following asymptotics: 



b(p) 



p for 1 — p <C 1, 
forp^p trG0 . 



(1) 



Clearly, a grey bond may never be removed, and at 
P = Ptrcc all bonds are either removed or grey. It is easy 
to see that the backbone of the entire (black and grey) 
network coincides with the backbone of the black sub- 
system. Indeed, no grey bond can belong to backbone: 
since its removal produces an isolated finite cluster, such 
a bond belongs to a dangling end. Note that black bonds 
are removed totally at random, hence the behavior of the 
black subsystem is identical to that of the standard bond- 
percolation system. In particular, the backbone vanishes 
at the percolation point, where b — p p0 rc- It follows that 
there exists a critical concentration of bonds p c such that 
for p < p c the remaining bonds all belong to one infinite 
cluster (which has finite density), while this cluster has 
no backbone. The critical concentration is determined by 
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FIG. 2: Typical configurations of bonds in the solid and tree- 
like phases (we do not distinguish black and grey bonds here), 
together with corresponding block graphs. Black square — the 
infinite block (the backbone); large solid circles — finite non- 
trivial blocks; small solid circles — trivial blocks; small open 
circles — articulation points. 



the condition 



KPc 



Ppcrc? 



(2) 



where J3 p0 rc is the percolation threshold for the standard 
bond percolation problem on the same lattice. On the 
other hand, the number of grey bonds in the system is 
equal to the number of finite black clusters: g — n c \ (see 
Fig. [IJ. Thus the critical concentration p c can be ex- 
pressed solely through the characteristics of the standard 
percolation problem: 



Pc = Ppc 



(3) 



where n*j is the number of finite clusters (per one bond 
of the initial lattice) at the critical point. The latter is 
known for many lattices; in particular, for the square 
5)/4 (see jg) and p pC rc = 1/2, so 



lattice = 



Pc 



3\/3-3 



0.54904. 



(4) 



In scaling theory of percolation the relations between 
different critical exponents are normally derived from 
considerations involving distribution function of finite 
clusters (see, e.g., |2j{). In our model finite clusters do 
not exist whatsoever, but fortunately one can introduce 
blocks — alternative objects, which to some extent play 
the role of finite clusters and make it possible to develop 
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the scaling theory. By definition, a block is a maximal 
subgraph that cannot be disconnected by deletion of a 
single vertex ||. It is not difficult to show that either 
a block consists of a single bond (and its two ends), or 
any two bonds belonging to a block lie on a common cy- 
cle. Two distinct blocks may have at most one point in 
common; such a point is called an articulation point, and 
its deletion necessarily disconnects the system. Given a 
network, one can form a graph with blocks and articula- 
tion points of the network as vertices, with two vertices 
connected if they correspond to an articulation point and 
a block that contains it. Such a block graph is always a 
tree; an example is shown in Fig. [21 

The backbone that exists in the solid phase constitutes 
the only infinite block in the system; it has finite density 
and contains infinite cycles. The backbone is linked to 
an infinite number of branches — dangling ends, each dan- 
gling end being a finite tree of finite blocks. 

In the tree-like phase the infinite block collapses, so 
that there are only finite cycles in this phase. The corre- 
sponding bond configurations we will call quasitrees. 

In the solid phase, the backbone becomes more loose 
as p approaches p c . The fraction of bonds belonging to 
it tends to zero: 



P B (p) cx (p- Pc ) 0E 



(5) 



where (3b is the index of the backbone density for the 
standard percolation (in particular, (3b ~ 0.48 for d = 2, 
see, e.g., @ and 

The total number of dangling ends decreases as p — * p c 
from above, while the number of blocks in each dangling 
end and the number of bonds in a typical block increase 
and diverge as p — > p c . It is convenient to introduce the 
distribution function of finite blocks consisting of I bonds: 



n;(p)~r T / \l{p-p c ) 1/a 



(6) 



where f(x) is a universal function that decays exponen- 
tially at x 3> 1. Application of standard scaling argu- 
ments 2] to blocks lead to the following relations between 
the critical exponents: 
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FIG. 3: The path, connecting points A and B on the block 
graph; nontrivial blocks are shown as grey islands. The inter- 
nal structure of a typical block is shown in inserts: a typical 
MST-path, traversing a block (lower panel); the shortest path 
across the block (upper panel). 



Finally, 



Ih 



T-2 



(9) 



where ds is the fractal dimension of the backbone at 

P = Pc- 

Since the conduction process involves only the back- 
bone, the conductivity of the SRBP model is identical to 
that of the standard percolation, 



csrbp(p) = o- perc [b{p)] oc (p-PcT, 



(10) 



hence the corresponding critical exponent /i is the stan- 
dard one. For the minimal path length in the solid phase 
one has 



£(R) = 



R 
v(p)' 



(11) 



V = 



ad 



f(p) ~ {p-PcY 



(7) 



The length £(p) characterizes correlations within the 
backbone, and since the backbone for the SRBP model 
is the same as for standard percolation, we conclude that 
the exponent v for the SRBP model coincides with that 
for standard percolation. 

The exponent 7 that characterizes the behavior of the 
mean size S(p) of finite blocks near p c is 



S{p) 



P-Pb(p) 



« (p-PcY 



(8) 



where d^??* ^ is the graph dimension for the infinite clus- 
ter at the critical point for the standard percolation 
problem. As usual, the formula 1)11(1 is valid only for 
R in the opposite case R <C £(p) it should be 

substituted by the critical law 



t{R) oc R d 



(12) 



For R ~ the expressions <|II|) and (|I2|I match. 

At p — ptrcc our system is reduced to the MST en- 
semble. In high dimensions d > d c a minimal spanning 
tree on an infinite lattice may in fact have many com- 
ponents. It is believed (see |f0|) that d c = 8, and in 
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FIG. 4: Root mean square euclidean displacement vs. chem- 
ical distance for different values of p, obtained via averag- 
ing over 1600 realizations of a 2048 x 2048 lattice with peri- 
odic boundary conditions. Asymptotic slopes correspond to 
dmin(p c ) = 1.13(1) and d min (p c -0.01 = 0.54) = d m in(ptree) = 
1.22(1) in agreement with known results. 



d < 8 dimensions for almost all trees of the MST ensem- 
ble there exists a unique finite path V(A, B) connecting 
any two given sites A and B. This path is a certain 
non self-intersecting random walk with fractal dimension 
df = d m i n . While this dimension is known exactly for the 
2D uniform spanning tree (UST) ensemble, only numeri- 
cal estimates are available for the MST case: <i m i n w 1.22 

a 1.59 for 



for d = 2, _4 



1.42 for d — 3, and d n 



d = 4 (see Below we demonstrate that the graph di- 

mension is the same throughout the entire tree-like phase: 



V(A,B)(Q) leading from A to B is, obviously, the mini- 
mal path over the graph V(A, B)^ . On nontrivial (con- 
taining more than one bond) blocks the path V{A, B)^ 
is the shortest path that crosses the block; it may be con- 
siderably shorter than any individual MST-path. This 
consideration enables one to estimate the typical ratio of 
lengths for a piece of the minimal path V(A, B)^ and 
the corresponding piece of the MST minimal path. We 
make such an estimate for the case when p < p c but 
p c — P <C 1 (i.e., for the vicinity of the phase transition). 
Having in mind that the typical block size is £(p) 3> 1, 
for a typical length MST-path crossing such a block we 

,(MST) 

get £mst{0 ~ ( ° iin , while the shortest path travers- 
ing the block is the same as for the critical percolation: 
4hort(0 ~ . As a result, we arrive at the estimate 

(|14|l . which matches with 1)12(1 for the critical case R ~ £. 

Since the above consideration is not quite rigorous, we 
have also undertaken numerical evaluation of d m i n (p) in 
order to check the identity 1(13(1 . Simulations did not 
show any variation of d m ; n with p (see Fig. 0}. 

In conclusion, we have demonstrated, both numerically 
and analytically, that the Self-Repairing Bond Percola- 
tion model undergoes a topological phase transition at a 
certain concentration p c of remaining bonds. In the tree- 
like phase (for p < p c ) the network, although being fully 
connected, has no backbone and hence zero conductiv- 
ity. The corresponding graphs of bonds are "quasitrees" : 
they contain only finite cycles (even for the infinite lat- 
tice) . The properties of the statistical ensemble of qua- 
sitrees are similar to those of the Minimal Spanning Trees 
ensemble. 



rfmin(p) = d\ 



(MST) 
min 



(13) 



More precisely, 



e(R, P ) 



R 



,(MST) 



v(p) 



v(p) ~ 



.(MST) ,(pcrc) 



,(MST) . ,(pcrc) 



(14) 

Consider a certain quasitree Q and the set of trees 
T} Q) which can originate from Q in the course of further 
destruction of bonds. Obviously, Q is a union of all : 

Q = urf Q) . (15) 
Now we introduce a graph 

V(A,B)^ = U i P i (A i B)W, (16) 

which is the union of all paths leading from A to B 
in all trees 7^ (Q) . It is easy to show that V(A,B) (Q ^ is 
precisely the path leading from A to B in the block graph 
(see Fig. |3Jl. The minimal (over the entire quasitree) path 
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